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14 Abstract. 

15 The 14-yr (1991-2005) time series of mesospheric water vapor from the Halogen Occultation 

16 Experiment (HALOE) are analyzed using multiple linear regression (MLR) techniques for their 

17 seasonal and longer-period terms from 45 S to 45N. The distribution of annual average water 

18 vapor shows a decrease from a maximum of 6.5 ppmv at 0.2 hPa to about 3.2 ppmv at 0.01 hPa, 

19 in accord with the effects of the photolysis of water vapor due to the Lyman-a flux. The 

20 distribution of the semi-annual cycle amplitudes is nearly hemispherically symmetric at the low 

21 latitudes, while that of the annual cycles show larger amplitudes in the northern hemisphere. The 

22 diagnosed 1 1-yr, or solar cycle, max minus min, water vapor values are of the order of several 

23 percent at 0.2 hPa to about 23% at 0.01 hPa. The solar cycle terms have larger values in the 

24 northern than in the southern hemisphere, particularly in the middle mesosphere, and the 

25 associated linear trend terms are anomalously large in the same region. Those anomalies are 

26 due, at least in part, to the fact that the amplitudes of the seasonal cycles were varying at northern 

27 mid latitudes during 1991-2005, while the corresponding seasonal tenns of the MLR model do 

28 not allow for that possibility. Although the 1 1-yr variation in water vapor is essentially 

29 hemispherically-symmetric and anti-phased with the solar cycle flux near 0.01 hPa, the 

30 concurrent temperature variations produce slightly colder conditions at the northern high 

31 latitudes at solar minimum. It is concluded that this temperature difference is most likely the 

32 reason for the greater occurrence of polar mesospheric clouds at the northern versus the southern 

33 high latitudes at solar minimum during the HALOE time period. 

34 

35 
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36 1 . Background and Objectives 

37 Water vapor (H 2 O) in the low to middle mesosphere is detennined by its so-called “entry level” 

38 value to the stratosphere from the tropical upper troposphere plus the effects of the oxidation of 

39 methane (CH 4 ) to H 2 O in the upper stratosphere. That conversion of CH 4 to H 2 O is essentially 

40 complete for air reaching the middle mesosphere. Water vapor and its component radicals also 

41 interact chemically with odd oxygen in the mesosphere. The distribution of H 2 O in the middle 

42 and upper mesosphere is determined by the effects of global-scale, net transport processes and by 

43 photolysis from Lyman-alpha (Ly-a) radiation that is most efficient near the mesopause. The 

44 largest periodic variations of TbO in the upper mesosphere are found at the higher latitudes and 

45 on both seasonal and 1 1-yr (or solar cycle) time scales. Thus, the observed variations of H 2 O 

46 can be a useful diagnostic of the perfonnance of radiative-chemical-transport models and of 

47 global change in the mesosphere [ Brasseur and Solomon, 1984; Lastovicka et a/., 2008]. 

48 

49 Currently, only a few multiyear datasets are available for characterizing the large-scale variations 

50 of water vapor in the mesosphere. In particular, H 2 O profiles were obtained with a vertical 

51 resolution of about 2.5 km from cloud tops to just above the 0.01-hPa level for more than one 

52 complete solar cycle (1991-2005) by the Halogen Occultation Experiment (HALOE) of the 

53 Upper Atmosphere Research Satellite (UARS) [ Russell et al., 1993; Harries et a/., 1996; Kiev et 

54 al., 2000]. Individual retrieved H 2 O mixing ratio profiles appear rather noisy (or oscillatory with 

55 altitude) in the mesosphere, due to the non-linear relationship between the retrieved H 2 O and the 

56 small changes of the observed 6.6-pm limb transmission values. Therefore, this study analyzes 

57 time series of points obtained by averaging mixing ratio profiles for each of the HALOE sunset 
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58 (SS) and sunrise (SR) orbital crossing occurrences of selected latitude zones. In order to obtain 

59 reasonable zonal mean estimates for a latitude crossing, the profiles are binned according to 20 

60 degree-wide latitude zones. A minimum of 5 profiles is required for a bin-average, although 

61 most of the averages are derived from many more profiles. These criteria are providing time 

62 series points that are much less noisy and are representative of zonal means because the large- 

63 scale longitudinal gradients of H 2 O are weak in the mesosphere at low and middle latitudes. 

64 

65 Section 2 describes briefly the multiple linear regression (MLR) model analysis approach that is 

66 used to obtain the seasonal and longer-term variations of HALOE Version 19 (VI 9) H 2 O versus 

67 latitude and for pressure-altitudes from 0.2 to 0.01 hPa. It also provides the latitude versus 

68 pressure distributions of the amplitudes of the semi-annual and annual tenns from the regression 

69 models. Section 3 shows results for the 1 1-yr or the solar cycle (SC-like) terms and relates the 

70 distribution of their “max minus min” H 2 O values with those from model simulations reported in 

71 the literature. It also contains findings for the concurrently analyzed, linear trend terms and of 

72 the potentially confounding effects for the detennination of their SC-like responses. Section 4 

73 discusses the SC-like variations in temperature and in TLO plus their relation to the occurrences 

74 of polar mesospheric clouds (PMC) near 0.01 hPa. Section 5 summarizes the findings. 

75 

76 2. Analysis Procedure and the Seasonal and QBO-like Variations 

77 Rems berg [2008a, his Figure 1] showed that the time series of tangent point measurement 

78 locations for the HALOE SR and SS events of 2001, for example, are representative of each of 

79 the seasons and for both hemispheres, at least Equatorward of about 55 degrees of latitude. 
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Conversely, the very high latitudes were not sampled at all during winter, due to the orbital 
inclination of UARS and of the Earth with respect to the Sun. Mesospheric H 2 O attains its 
minimum values during winter at high latitudes, which means that the amplitude of its annual 
cycle may be significantly underestimated from the HALOE data in those regions of the upper 
mesosphere. It was also noted by McHugh et al. [2003] that a PMC feature in the measurement 
line-of-sight can contaminate the retrieved H 2 O profile significantly. Time series of HALOE 
H 2 O data are not considered for analysis at the higher latitudes for these reasons. 


The MLR analysis approach of the present study follows that given in Remsberg [2008a] and is 
described briefly here. Zonal, bin-averages of the V19 SS and then the SR profiles are obtained 
for ten, 20-degree wide latitude zones from 45S to 45N and for 9 pressure levels from 0.2 hPa (or 
near 60 km) to 0.01 hPa (or near 80 km). Each of the 90 time series of the set of SS plus SR H 2 O 
is analyzed by MLR techniques to obtain their annual oscillation (AO) and semi-annual 
oscillation (SAO) terms, any weak quasi-biennial oscillation (QBO)-like and sub-biennial 
(interannual or IA) tenns, their decadal-scale (1 1-yr) terms, and linear trend terms. The sub- 
biennial term occurs as a result of interactions between the AO and QBO tenns, and it has a 
period of about 21 months. For a given latitude bin and pressure-altitude the amplitudes and 
phases of the periodic terms, including the QBO and 1 1-yr terms, are based on a fitting with a set 
of predictors that are simply hannonics of specific periods. This approach differs from the more 
common practice of fitting to a proxy for an underlying physical process, such as the fl0.7-cm 
flux for the SC and a couple of tropical wind indices for the QBO tenn (e.g., Randel et al, 2000; 
Wallace et al., 1993). On the other hand, the present approach is particularly useful for 
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102 identifying interannual responses in the mesosphere and for detennining whether the primary 

103 decadal-scale response of H 2 0 is truly anti-phased with the SC uv-flux forcing. 

104 

105 It is easy to confound the effects of SC and trend terms for the time series, especially when 

106 adjacent zonal mean points are correlated. Those effects are accounted for to first order by 

107 conducting a two-step process with the MLR model in the manner of Tiao et al. [1990] and 

108 Rems berg [2008a], Initially, the set of relevant model terms are fit to the FLO time series, and a 

109 weak positive, lag-1 autoregression (AR1) coefficient is obtained from the time series of the 

110 noise residuals. Then, the model tenns are transfonned to account for the AR1 coefficient, and 

111 the data time series are refit to obtain the coefficients for the final MLR terms. 

112 

113 At this point it is noted that the 14-yr time series consist of alternating SR and SS points and that 

114 there are slight differences between the means of the SR and SS H 2 O values for several extended 

115 periods of the dataset. Those differences are due to small biases from the FOV lockdown of the 

116 HALOE instrument in its SR transmission profiles. The effect is evident as an exaggerated 

117 vertical oscillation in the retrieved SR H 2 O profile near 0. 1 hPa — the FLO values are too small 

118 just below and too large just above that pressure level. The incorrect lockdown for those SR 

119 scans leads to weakly negative AR1 coefficients that are accommodated by the transfonned 

120 tenns of the final MLR model for the combined (SR plus SS) HALOE FLO time series. 

121 
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122 Figure 1 is an example time series of HALOE VI 9 FF>0 at 35N and 0.015 hPa (near 75 km). The 

123 solid and open points are the bin-averaged SS and SR values, respectively, for HALOE tangent- 

124 point measurements made at 35±10N throughout the years. The oscillating curve is the MLR 

125 model fit to the points, and it is composed of constant (Const), AO, SAO, QBO (853-dy or 28 

126 mo), IA (640-dy or 21 -mo), 1 1-yr (or SC-like), and linear (Lin) trend terms. There is a clear 

127 annual cycle in the data at this latitude and level. The straight line in Figure 1 is just the sum of 

128 the Const and Lin tenns. Figure 2 is the time series of the data minus model residuals (in ppmv) 

129 from Figure 1. No apparent periodic structure is remaining in those residuals — an important test 

130 for the acceptance of the final MLR model. 

131 

132 The HALOE HzO data time series of Figure 1 exhibits a weak, near-biennial variation; the 

133 amplitude of the associated, model QBO (28-mo) tenn is discussed later in this section. There 

134 are also clear relative minima at about 1991 and 2002, in accord with the estimated times for the 

135 maxima of the Ly-a photolysis flux of the 1 1-yr solar cycle. Such a direct response is expected 

136 for the upper mesosphere. The response to Ly-a is much weaker in the middle mesosphere and 

137 is most likely mixed with the larger H 2 O responses of the higher altitudes and latitudes as a result 

138 of its annually-averaged, net circulation. In order to assess the prospect of a delayed, SC-like 

139 response for the H 2 O of the middle mesosphere, a simple sinusoid of 1 1-yr period was used to fit 

140 the data at all altitudes and latitudes. The amplitudes of the 1 1-yr terms were determined, but 

141 then the phases were checked to see how well the 1 1-yr H 2 O maxima coincided with the 

142 estimated minima for a Ly-a or other solar cycle (SC) flux proxy. Findings for the 1 1-yr terms 

143 are presented and discussed in Section 3. 
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Figures 3-5 show the latitude versus pressure cross sections of (a) the constant term (in ppmv), 

(b) the amplitude of the SAO tenn (in %), and (c) the amplitude of the AO tenn (in %), 
respectively. Since HALOE samples a latitude zone infrequently, one must combine the 
constant tenn, seasonal amplitudes, and their phases, in order to generate an H 2 O distribution for 
comparisons with models at a specific time of the year, for example in January. Jackson et al. 
[1998] showed and discussed these tenns for the mesosphere from the first 5 years of the 
HALOE data, and the results of this section are compared with them briefly. Figure 3 is the 
distribution of the constant tenns, i.e., their 14-yr, annual average H 2 O mixing ratios. Those 
values decrease from a maximum of about 6.5 ppmv in the mid-mesosphere to 3.2 ppmv near 
0.01 hPa. Estimates of total bias enor for HALOE zonal mean H 2 O are of order 10% in the 
mesosphere [Harries et al., 1996]. Comparisons with other satellite data indicate that the 
HALOE values may be, in fact, too low but by no more than 12% in the same altitude region 
[Lambert et al, 2007; Milz et al., 2009]. The relative variations of HALOE H 2 O with altitude 
and latitude are likely more accurate because the primary systematic errors do not depend on 
atmospheric state. The decreasing values of H 2 O with altitude are due largely to the effects of 
Ly-a photolysis. Their decrease toward higher latitudes is indicative of the descent of air at polar 
winter latitudes. There is also a slight hemispheric asymmetry for the northern versus the 
southern latitudes for this distribution. 


Figure 4 is the distribution of the SAO amplitudes, which have their greatest percentage values 
(or order 20%) near 0.01 hPa. This distribution agrees well with the observations of the Sub- 
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166 Millimetre Radiometer (SMR) instrument of the ODIN satellite [ Lossow et al., 2008], even near 

167 0.01 hPa, where the HALOE fEO profiles are close to their signal-to-noise limit. The HALOE 

168 SAO amplitudes occur nearly symmetric about the Equator and with their minimum values in the 

169 subtropics. The distribution of the phases of the SAO maxima for their first cycles is shown in 

170 Table 1 (in days from January 1) and is also nearly hemispherically-symmetric. A steady 

171 descent is indicated for the SAO from the upper to the middle mesosphere at the low latitudes 

172 (15S to 15N), in accord with the SAO phase observations and forcings for other atmospheric 

173 parameters [ Garcia et al., 1997]. At the middle latitudes the descent of the SAO is also 

174 apparent in the southern hemisphere, while it is variable or nearly stationary with altitude in the 

175 northern hemisphere. 

176 

177 Figure 5 is the distribution of the amplitudes of the annual cycle (AO) terms; amplitudes exceed 

178 30% at the high altitudes and latitudes. The hemispheric asymmetry of this tenn is quite striking 

179 and is responsible for the observed hemispheric asymmetry in the annual average mixing ratios 

180 of Figure 3. Minimum amplitudes occur along a vertical axis centered near 15S, and this finding 

181 agrees very well with the SMR observations of Lossow et al. [2008, see their Figure 8]. Similar 

182 plots in Jackson et al. [1998] show that this vertical axis is nearer to 5S. On the other hand, the 

183 AO amplitudes in Figure 5 are approximately the same at 45S and 45N at the 0.2-hPa level. The 

184 generally smaller AO amplitudes of the southern hemisphere may be due, in small part, to the 

185 fact that the Earth-Sun distance is at a minimum and the associated Ey-a flux at a maximum in 

186 SH summer, which would slightly offset the tendency of the summer upwelling to increase TEO. 

187 However, this circumstance cannot explain the axis of minimum amplitudes at 15S. 
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188 


189 Table 2 shows the time of occurrence of the AO maximum in FLO in days past January 1 . Those 

190 phases are anti-symmetric between the two hemispheres (but seasonally symmetric), especially 

191 at the middle latitudes. Phase propagation is more complicated in the subtropics, however. For 

192 example, at the 0.05-hPa level the AO phase at 15N is very similar to that of northern middle 

193 latitudes, while the AO phase at 15S is one-quarter cycle out-of-phase with that at 35-45S. 

194 Clearly, the seasonal cycle in the NH versus the SH is projecting differently onto the 12-mo and 

195 6-mo hannonics of the MLR fit, especially for the subtropics. Those two terms have significant 

196 amplitudes at 15N and 0.05 hPa as indicated in Figure 6, while the structure of the points in the 

197 corresponding time series at 15S is much weaker (not shown). Most likely the characters of 

198 both the AO and the SAO in the middle mesosphere are affected by planetary Rossby waves 

199 l Dunkerton and Delisi , 1985] and/or gravity waves [Sato et a/., 2009], which are 

200 hemispherically-asymmetric in their occurrence. 

201 

202 Figure 7 is the distribution of the amplitudes of the QBO-like terms. Largest amplitudes occur at 

203 middle latitudes at about 0.03 hPa, and another maximum occurs in the tropical upper 

204 mesosphere. In fact, the distribution in Figure 7 is very similar to that of the QBO terms for the 

205 HALOE temperature data in Remsberg [2008b], indicating that the interannual variability of FLO 

206 is primarily due to transport processes. The QBO amplitudes are somewhat larger in the 

207 northern than in the southern hemisphere subtropics, and the minimum values occur at about 

208 15S. These differences may indicate that the wave forcings associated with at least one phase of 
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209 the QBO cycle are interacting with the SAO forcings of the northern subtropics more effectively 

210 [Garcia et al., 1997]. 

211 

212 At this point it is noted that the MLR analyses herein were conducted independently for each of 

213 the 90 HALOE H 2 O time series. Therefore, an important measure of the quality of the results in 

214 Figures 3-5 and 7 is the spatial coherence of their patterns with altitude and latitude. Uncertainty 

215 for the coefficients of the tenns from the MLR fit to the data is of the order of 0.20 ppmv (2-o), 

216 which is small compared with the amplitudes of the seasonal terms. Distributions for the QBO 

217 and IA terms are not very significant because their amplitudes are less than 0.22 ppmv (4%) 

218 everywhere. One should also be aware that the estimates of uncertainty assume that the various 

219 tenns of the MLR models are orthogonal to each other. That requirement is not met for the 1 1-yr 

220 and linear trend terms of the 14-yr HALOE time series. 

221 

222 3. Decadal-scale Responses in Water Vapor 

223 Early analyses of the mesospheric response of HALOE water vapor to the variations of solar Ly- 

224 a were provided by Randel et al. [2000] and Chandra et al. [1997]. Randel et al. [2000] 

225 reported a “max minus min” variation of 2 1% for the period 1992-99, while Chandra et al. 

226 [1997] reported an increase of 41% from January 1992 to December 1995, both at 0.01 hPa. Part 

227 of their differences may be traceable to the HALOE version 19 dataset used by Randel et al. 

228 [2000] versus the version 18 data employed by Chandra et al. [1997]. Their findings may also 

229 have been affected by the underlying trends for H 2 O that are increasing in the mesosphere 

230 through the 1990s; the trends are more nearly flat thereafter [ Nedoluha et al., 2009]. 
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231 


232 As noted in Section 2, the approach to the SC analyses herein has been to simply fit the HALOE 

233 water vapor time series with a sinusoid of 1 1-yr period and then to check its phase to see how 

234 closely it conforms to an estimate of the variations of the direct Ly-a flux. Since EEO is 

235 photo lyzed most effectively at solar maximum, its 1 1-yr responses are plotted for solar “max 

236 minus min” (as a percent of its average value) and are therefore negative. Figure 8 is a contour 

237 plot of that full change in H 2 O (i.e., twice its amplitude), which varies from about 4% near 0.2 

238 hPa to about 23% at 0.01 hPa, the latter value agreeing most closely with the values of Randel et 

239 al. [2000] and Nedoluha et al. [2009]. Values of 4% are just significant at the 2-a level; the 

240 greater percentage values are considered highly significant. 

241 

242 Figure 9 is the phase of the maximum of the 1 1-yr PFO response (in years past January 1991 or 

243 2002). The dashed contour denotes the 5.5-yr value or July 1996, the approximate time of solar 

244 flux minimum (or water vapor maximum). In the middle mesosphere the water vapor maximum 

245 is lagging solar flux minimum by about 1 year at middle latitudes and by nearly 2 years in the 

246 tropics. However, the phases of those weak amplitude terms are just barely significant, and they 

247 are the tenns that are most likely to be confounded with the associated linear trend terms (see 

248 further discussion in the last paragraph of this section). 

249 

250 Figure 10 is the plot of the so-called SC-like HoO responses. They were obtained by adjusting 

251 the 1 1-yr responses of Figure 8 by the fact that they were not quite in-phase with the Fy-a flux 
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252 forcing — i.e., an adjustment by the cos [(t-5.5)/l 1], where t is the time of the 1 1-yr water vapor 

253 maximum from Figure 9. In effect, this adjustment accounts for the possibility that there are 

254 other decadal-scale forcings for the H 2 O that were not considered in the MLR model. There is 

255 little change for the H 2 O responses between Figures 10 and 8 for the upper mesosphere. Those 

256 max minus min FLO percentage variations are opposite in sign but very similar in magnitude to 

257 the observed 1 1-yr changes in the Lyman-a flux [ Tobiska et al., 1997; De/and et al., 2004]. The 

258 adjusted, max minus min FLO variations are altered the most for the tropical middle mesosphere 

259 and are now no larger than about 2%. The 1 1-yr and the SC-like responses are nearly symmetric 

260 about the Equator in the uppennost mesosphere. However, the diagnosed responses at middle 

261 latitudes are about twice as large in the northern as in the southern hemisphere, indicating a 

262 potentially confounding influence from the trend tenns. 

263 

264 Model H 2 O solar cycle responses agree well with the HALOE findings in the uppermost 

265 mesosphere, once an accommodation is made for the fact that the percentage responses for many 

266 of the models are referenced to their 1 1-yr H 2 O minimum values rather than to the response 

267 average, as in Figures 8 and/or 10 [e.g., Garcia et al., 1984; Huang and Brasseur, 1993; and 

268 Schmidt et al., 2006]. The SC-like response profile from Figure 10 at 20N also agrees well with 

269 that from the Water Vapor Millimeter-wave Spectrometer (WVMS) instrument at Mauna Loa for 

270 the concurrent period of the HALOE data [Nedoluha et al., 2009]. 

271 
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272 The MLR models include linear trend terms, and Figure 1 1 is their distribution (in %/decade). 

273 Values vary from 12 (or 1.2%/yr) at 40N and 8 (or 0.8%/yr) at 40S to nearly no change at the 

274 Equator. Such large values for the middle latitudes of the mesosphere are unphysical compared 

275 with the model simulations of about 4%/decade, due primarily to the trends from the oxidized 

276 methane [ Garcia et al, 2007]. It is also noted that no significant changes in the HALOE 

277 instrument have been found that might be contributing to these observed FLO trends and, in 

278 particular, their variations with latitude in Figure 1 1 [ Gordley et al., 2009]. 

279 

280 The distribution of trends in the tropics to subtropics of Figure 1 1 is qualitatively consistent with 

281 the FLO and CFL trends of the stratosphere of several years earlier. For example, Randel et al. 

282 [2006] reported FLO increases of order 1%/yr in the lower stratosphere from HALOE data for 

283 the decade of the 1990s, followed by a decreasing to no trend after 2001. Scherer et al. [2008] 

284 found HiO trends from balloon-borne, frostpoint measurements at Boulder of -0.2 to 1.0%/yr for 

285 altitudes from 14 to 25 km for the same period with a clear change of trend in 2000/2001. 

286 Trends in tropospheric CH 4 were variable and generally positive (~0.3%/yr) through the 1990s, 

287 but also slowed to near zero thereafter [Dlugokencky et al., 2009]. Correspondingly, the 

288 observed positive trends of CH 4 for the lower stratosphere slowed to about zero during the 

289 HALOE period [Rohs et al., 2006], It is reasonable to assume from the age-of-air characteristics 

290 for the H2O and CH4 entering the lower tropical stratosphere that their trends will be retained (at 

291 least to some extent) as that air is transported upward over a period of several years to near the 

292 stratopause by the Brewer/Dobson circulation [ Shepherd , 2007]. However, because the 

293 exchange of air between low and high latitudes is much more efficient in the mesosphere, the 
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294 diagnosed HALOE trends of the middle latitudes cannot be due to the changing stratospheric 

295 source gases. 

296 

297 The character of the MLR fit at northern middle latitudes is examined further in Figure 12 for 

298 35N and the 0.05-hPa level. A closer look at the annual cycle variations in the data reveals that 

299 they are smaller than those from the MLR fit during 1995-97, while they are clearly larger in the 

300 data from about 2000-03. In other words, the amplitudes of the seasonal cycles in the ELO data 

301 are varying over the 14-yr time period of HALOE. Thus, the residuals that are obtained after the 

302 time series data are de-seasonalized will contain some longer-period structures that are being 

303 accounted for by the interannual, 1 1-yr, and linear trend tenns of the MLR model. The 

304 constraint of the seasonal-cycle fit is contributing to the anomalous character of the diagnosed 

305 1 1-yr and trend terms. 

306 

307 It is postulated that the rather marked difference for the 1 1-yr and trend terms of the two 

308 hemispheres is also an indicator of the differences for the seasonal exchanges of air between the 

309 middle and high latitudes. During winter there is descent of air having low H 2 O values in the 

310 polar vortex. The degree to which that air is transported and/or mixed to lower latitudes depends 

311 on the extent to which the vortex remains intact through the winter/spring period, and generally 

312 there is more wintertime wave activity propagating to the mesosphere from below in the northern 

313 hemisphere [Sato et al., 2009]. Northern hemisphere, mid-winter warming events in the upper 

314 stratosphere are accompanied by a cooling in the mesosphere, a weakening or breakdown of the 
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315 polar night jet, and an enhanced meridional transport and mixing of the air. Such mid-winter 

316 events are rarely found in the southern hemisphere. Pawson and Naujokat [1999] reported that 

317 there was very little mid-winter wave activity at NH high latitudes during most of the 1990s, 

318 while Manney et al. [2005] found such activity to be prevalent for 6 of the next 7 years. Since 

319 HiO has a significant meridional gradient in the middle mesosphere at the mid latitudes (Figure 

320 1), it is reasonable to conclude that the degree to which the seasonal cycles of the model do not 

321 match the data in Figure 12 is an important indicator of the response of the transport to those 

322 events. Such activity may bear little relation to the solar cycle forcing, although there is no way 

323 to know for sure based on this rather short, 14-yr HALOE dataset. 

324 

325 Observations of increasing lower stratospheric H 2 O trends up to 2001, followed by decreasing 

326 trends thereafter, strongly suggest that 14 years is much too short a time for an analysis of its 

327 trends. The findings herein from the HALOE dataset also support that contention for the 

328 mesosphere. As Garcia et al. [2007] argued, one may need 30-40 years of data before one can 

329 estimate a true secular trend in water vapor for the middle atmosphere. 

330 

331 4. Implications of the Solar Cycle Responses in H 2 0 

332 

333 There have been a number of observational and model studies of the effects of water vapor in the 

334 upper mesosphere and its relation to the occurrence of PMCs near 0.01 hPa [e.g., Shettle et al, 

335 2009; Sonnemann and Grygalashvyly, 2005]. Although the analyses of the HALOE water vapor 

336 responses herein do not quite extend to the high latitudes of the PMCs, Figure 5 indicates that the 
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337 AO term has larger percentage responses in the northern than in the southern hemisphere for that 

338 region of the upper mesosphere. Maximum H 2 O mixing ratios occur in the summer at high 

339 latitudes. Recall that Figure 9 indicates that PMCs should also occur most frequently at solar 

340 minimum, when H 2 O values are at their maximum. On the other hand. Figure 10 shows that the 

341 percentage changes for the SC-like responses of H 2 O are nearly identical for the two 

342 hemispheres near 0.01 hPa, indicating that there should be no hemispheric difference in the SC- 

343 like occurrence of PMCs due to water vapor. 

344 

345 The observed greater seasonal occurrence of PMCs in the northern versus the southern 

346 hemisphere would be aided by slightly colder temperatures near the summer polar mesopause of 

347 the northern hemisphere [Wrotny and Russell, 2006; Hervig and Siskind; 2006]. Figure 13 is the 

348 plot of the HALOE AO temperature amplitudes for 10-degree wide latitude bins from 60S to 

349 60N adapted from Rems berg [2007]. Note that although the seasonal sampling from HALOE is 

350 barely adequate at the higher latitudes, signal-to-noise values in the mesosphere are much larger 

351 for the temperature than for the H 2 O measurements. Therefore, the MLR analyses were 

352 extended to higher latitudes for temperature, and it was found that its AO amplitudes agree 

353 closely with those from the earlier satellite temperature climatology of Barnett et al. [1985]. 

354 

355 Figure 13 shows that the seasonal temperature variations of the upper mesosphere extend farther 

356 into the subtropics for the northern hemisphere — an indication of the associated effects of the net 

357 transport for the two hemispheres. Note also that the AO temperature amplitudes at the high 
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358 latitudes are larger by up to 4 K in the southern than in the northern hemisphere, potentially 

359 leading to formation of more PMC in southern hemisphere summer. But when those AO 

360 temperature variations are referenced to the annually-averaged temperatures of Figure 14 that are 

361 colder by up to 5 K near 0.01 hPa in the northern versus the southern high latitudes, it is the high 

362 latitudes of the northern hemisphere that are actually favored slightly for PMC fonnation. Thus, 

363 it should be clear that if one wants to reconstruct the actual temperature from the MLR fit, one 

364 needs to know the constant tenn and resolve both the annual and semiannual harmonics from the 

365 HALOE data time series. Although the seasonal sampling for the HALOE FLO measurements 

366 was marginal at the high latitudes where PMCs occur, the combination of slightly colder 

367 temperatures plus slightly larger FLO mixing ratios at the nearby latitudes favors the formation 

368 of PMC at the northern high latitudes in summer, in accord with earlier findings [ Wrotny and 

369 Russell, 2006; Hervig and Siskind, 2006]. 

370 

371 Remsberg [2008b] found that the HALOE temperatures are also colder at solar minimum than at 

372 solar maximum at the higher latitudes near 0.01 hPa. In other words, when the SC-like H 2 O 

373 mixing ratios are at their maximum, the temperatures are coldest. Furthennore, the SC-like, max 

374 minus min, temperature responses (though small) are about a factor of 2 greater at the northern 

375 than at the southern high latitudes of the upper mesosphere, as shown in Figure 15. Therefore, it 

376 is tentatively concluded that the occurrence of PMCs is significantly greater at solar minimum in 

377 the northern hemisphere than in the southern hemisphere due to the temperature differences, in 

378 accord with the earlier estimates of Hervig and Siskind [2006], Siskind et al. [2005] and Wrotny 

379 and Russell [2006], Even so, it may be that this finding is unique to the HALOE time period, as 
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380 opposed to being a regular occurrence [e.g., Luebken et al., 2009]; longer time series are needed 

381 to verify that prospect. 

382 

383 Marsh et al. [2003] analyzed HALOE ozone for its response to HOx (and fEO) in the upper 

384 mesosphere and reported that only its SS (not SR) ozone was sensitive to those changes. They 

385 also found significant decreases in HALOE ozone associated with a nearly linear increase of 

386 about 1%/yr in the HALOE H 2 O, at least for the period of 1991-2001. However, the results of 

387 Figure 1 1 indicate much weaker H 2 O trends at most latitudes based on the full 14-yr HALOE 

388 dataset. It is very likely that the associated ozone trends are smaller, too. The overall response 

389 of mesospheric ozone to the H 2 O from HALOE should be reanalyzed and compared with the 

390 model simulations of the SC-like responses of ozone [e.g., Schmidt et al., 2006; Marsh et al., 

391 2007; Tsutsui et al., 2009]. Still, the direct effects of the solar cycle forcing for ozone are much 

392 smaller than for H 2 O in the upper mesosphere. Thus, the attribution of trends in ozone to any 

393 long-term changes in the H 2 O is affected by the precision of the ozone and requires further study. 

394 

395 5. Summary Findings 

396 The HALOE measurements were obtained with a sampling frequency that is adequate for 

397 resolving the seasonal and longer-tenn variations at low and middle latitudes. Its 14-yr (1991- 

398 2005) time series of H 2 O were analyzed for those variations in the middle and upper mesosphere. 

399 The distribution of the annual average H 2 O shows a decrease from a maximum of about 6.5 

400 ppmv in the mid-mesosphere to 3.2 ppmv near 0.01 hPa, in accord with the effects of the 

401 photolysis of H 2 O due to the Ly-a flux. Although SAO cycle amplitudes are hemispherically- 
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402 symmetric at low latitudes, AO cycle amplitudes are clearly larger in the northern than the 

403 southern hemisphere. It is presumed that this hemispheric asymmetry is reflective of a stronger 

404 net circulation for the northern latitudes, particularly during winter and springtime. 

405 

406 The response of a periodic, 1 1-yr term is essentially anti-phased with that of the solar cycle flux 

407 in the upper mesosphere. Its “max minus min” TEO responses vary from about 4% at 0.2 hPa to 

408 about 23% at 0.01 hPa; they are nearly hemispherically-symmetric in the uppennost mesosphere. 

409 At 20N the profile of the solar cycle (or SC-like) H 2 O response from HALOE is in good 

410 agreement with the microwave measurements of Nedoluha et al. [2009] and with numerical 

411 model simulations generally. The 1 1-yr response in the mid-mesosphere is nearly anti-phased, 

412 lagging that of the solar cycle flux by no more than 1 to 2 years; the response is asymmetric 

413 between the northern and southern hemispheres, however. The associated linear trend terms in 

414 the mid-mesosphere agree with model estimates at low latitudes, but are too large at the middle 

415 latitudes of the northern hemisphere. It was found that the observed seasonal cycle in TEO had 

416 varying amplitude over that 14-yr period that was not accounted for by the AO terms of the MLR 

417 model. The accompanying, de-seasonalized residuals led to anomalous results for the 1 1-yr and 

418 trend tenns in the northern hemisphere, as a result. Those variations of the observed seasonal 

419 cycle amplitudes are considered an indicator of changes in the meridional transport and mixing 

420 due to the wave forcing from below, although it is unclear whether there is any association with 

421 the solar cycle itself. The distribution of SC-like and trend terms may be unique to the 1991- 

422 2005 timeframe; data over several complete solar cycles are needed in order to know for sure. 
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423 Taken together, these analyzed, seasonal and decadal-scale H 2 O variations may be useful 

424 diagnostics for assessing the validity of model transport in the mesosphere. 

425 

426 The combination of high H 2 O and low temperature supports a higher frequency of occurrence of 

427 PMCs at high latitudes of the upper mesosphere during the summer. Although the SC-like 

428 variations in the HALOE H 2 O are nearly hemispherically-symmetric in the upper mesosphere, 

429 the concurrent variations in temperature are not. Temperatures are coldest at the higher latitudes 

430 of the northern hemisphere in summer and at solar minimum. It is concluded that these 

431 conditions are likely the primary reason for a greater occurrence of PMCs in the northern versus 

432 the southern hemisphere high latitudes at solar minimum during the HALOE time period. 

433 
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Table 1 — Phase of Maximum of First SAO Cycle (in days since Jan 1) 


600 


P(hPa) 

45S 

35 

25 

15 

5 

5 

15 

25 

35 

45N 

0.01 

5 

8 

10 

15 

32 

21 

20 

20 

14 

8 

0.015 

6 

5 

6 

23 

38 

39 

39 

32 

11 

8 

0.02 

10 

6 

181 

60 

71 

69 

72 

66 

6 

6 

0.03 

13 

10 

180 

95 

98 

96 

95 

107 

171 

3 

0.05 

24 

20 

12 

111 

115 

114 

113 

144 

122 

1 

0.07 

29 

27 

20 

150 

127 

141 

119 

144 

165 

5 

0.10 

39 

33 

26 

171 

140 

148 

143 

111 

8 

13 

0.15 

47 

37 

14 

173 

160 

159 

147 

111 

27 

8 

0.20 

58 

39 

5 

172 

171 

11 

157 

172 

20 

11 
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603 Table 2 — Phase of Maximum of Annual Cycle (in days since Jan 1) 

604 


P(hPa) 

45S 

35 

25 

15 

5 
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15 

25 

35 

45N 

0.01 

3 

3 

362 

4 

145 

157 

170 

180 

186 

185 

0.015 

10 

10 

2 

27 

120 

142 

167 

186 

193 

194 

0.02 

11 

14 

11 

40 

117 

141 

173 
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0.03 

12 

16 

29 

54 
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157 

187 
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0.05 

16 

21 

49 
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213 
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211 

0.07 

17 

21 

34 

115 
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199 

217 
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0.10 

15 

23 

27 

80 

175 

183 

201 

221 

224 

224 

0.15 

20 

22 

10 

43 
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179 
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224 

0.20 

21 

21 

6 

362 

185 

166 
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206 

231 

231 
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607 Figure 1. Time series of bin-averaged HALOE SS (solid) and SR (open) H 2 O values at 35N and 

608 0.015 hPa (near 75 km). The oscillating curve is the MLR model fit to those values. 

609 

610 Figure 2. Time series of the data minus model H 2 O residuals (in ppmv) of Figure 1. 

611 

612 Figure 3. Pressure versus latitude contour plot of the 14-yr, annual average H 2 O mixing ratio (in 

613 ppmv). Contour interval is 0.4 ppmv and altitude coordinate is approximate. 

614 

615 Figure 4. As in Figure 3, but for the amplitudes of the semi-annual (SAO) terms as a percentage 

616 of the 14-yr average values. Contour interval is 5%. 

617 

618 Figure 5. As in Figure 3, but for the amplitudes of the annual (AO) terms as a percentage of the 

619 14-yr average values. Contour interval is 5%. 

620 

621 Figure 6. As in Figure 1, but for 15N and 0.05 hPa. 

622 

623 Figure 7. As in Figure 3, but for the amplitudes of the quasi-biennial (QBO) terms as a 

624 percentage of the 14-yr average values. Contour interval is 0.5%. 
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625 Figure 8. As in Figure 3, but for the 1 1-yr “max minus min” responses, as a percentage of the 

626 14-yr average values. Contour values are negative and their interval is 2 %. 

627 

628 Figure 9. Time of the maximum of the 1 1-yr term for FFO (in years past January 1991 or 2002). 

629 The dashed curve denotes 5.5 years (or July 1996), and the contour interval is 0.2 years. 

630 

631 Figure 10. As in Figure 8, but for the solar cycle (or SC-like) “max minus min” H 2 O responses 

632 as a percentage of the 14-yr average values. Contour interval is 2%. 

633 

634 Figure 11. Pressure versus latitude plot of the linear trends for H 2 O in terms of their percentage 

635 changes per decade. The zero contour is dashed, and the contour interval is 1.5%/decade. 

636 

637 Figure 12. As in Figure 1, but for 35N and 0.05 hPa. 

638 

639 Figure 13. Pressure versus latitude plot of the temperature amplitudes for the annual cycle (in 

640 K). Contour interval is 2 K. 

641 

642 Figure 14. As in Figure 13, but for the 14-yr, annual average temperatures (in K). Contour 

643 interval is 5 K. 
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644 


645 Figure 15. As in Figure 13, but for the solar cycle (or SC-like), max minus min, temperatures 

646 (in K). The negative and zero contours are dashed, and the interval is 0.5 K. 

647 
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648 

649 Figure 1. Time series of bin-averaged HALOE SS (solid) and SR (open) H 2 O values at 35N and 

650 0.015 hPa (near 75 km). The oscillating curve is the MLR model fit to those values. 
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651 

652 Figure 2. Time series of the data minus model H 2 O residuals (in ppmv) of Figure 1. 

653 
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654 

655 Figure 3. Pressure versus latitude contour plot of the 14-yr, annual average H 2 O mixing ratio (in 

656 ppmv). Contour interval is 0.4 ppmv and altitude coordinate is approximate. 
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659 

660 Figure 4. As in Figure 3, but for the amplitudes of the semi-annual (SAO) terms as a percentage 

661 of the 14-yr average values. Contour interval is 5%. 
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667 

668 Figure 6. As in Figure 1, but for 15N and 0.05 hPa. 
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670 

671 Figure 7. As in Figure 3, but for the amplitudes of the quasi-biennial (QBO) terms as a 

672 percentage of the 14-yr average values. Contour interval is 0.5%. 
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673 

674 Figure 8. As in Figure 3, but for the 1 1-yr “max minus min” responses, as a percentage of the 

675 14-yr average values. Contour values are negative and their interval is 2%. 
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681 

682 Figure 10. As in Figure 8, but for the solar cycle (or SC-like) “max minus min” H 2 O responses 

683 as a percentage of the 14-yr average values. Contour interval is 2%. 
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686 

687 Figure 11 . Pressure versus latitude plot of the linear trends for H 2 O in terms of their percentage 

688 changes per decade. The zero contour is dashed, and the contour interval is 1.5%/decade. 
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690 

691 Figure 12. As in Figure l, but for 35N and 0.05 hPa. 
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694 

695 Figure 13. Pressure versus latitude plot of the temperature amplitudes for the annual cycle (in 

696 K). Contour interval is 2 K. 
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699 Figure 14. As in Figure 13, but for the 14-yr, annual average temperatures (in K). Contour 

700 interval is 5 K. 
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702 

703 Figure 15. As in Figure 13, but for the solar cycle (or SC-like), max minus min, temperatures 

704 (in K). The negative and zero contours are dashed, and the interval is 0.5 K. 
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